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ABSTRACT 

We investigate the structure of dark matter halos by means of the kinematics of a very large sample 
of spiral galaxies of all luminosities. 

The observed rotation curves show an universal profile which is the sum of an exponential thin disk term 
and a spherical halo term with a flat density core. We find that the Burkert profile proposed to describe 
the dark matter halo density distribution of dwarf galaxies also provides an excellent mass model for the 
dark halos around disk systems up to 100 times more massive. Moreover, we find that spiral dark matter 
core densities po and core radii tq lie in the same scaling relation po — 4.5 x lQ~'^{r{i/kpc)~'^/^MQpc~^ 
of dwarf galaxies with core radii up to 10 times smaller. 

At the highest masses po decreases with rg faster than the — | power law implying a lack of objects 
with disk masses > IO^^Mq and central densities > 1.5 x 10^^(ro/fcpc)^'^ M^pcr^ that can be explained 
by the existence of a maximum mass of ~ 2 x lO^^M© for an halo hosting a spiral galaxy. 

Subject headings: Galaxies: spiral, kinematics and dynamics Cosmology: dark matter 



1. INTRODUCTION 

It is now well established that spiral galaxies have uni- 
versal rotation curves (URC) that can be characterized 
by one single free parameter, the luminosity (e.g. Ru- 
bin et al. 1980; Persic & Salucci, 1991, Persic, Salucci & 
Stel 1996 (PSS)). For instance, low-luminosity spirals show 
ever-rising rotation curves (RC) out to the optical radius 
^, while, in the same region, the RC of high-luminosity 
spirals are flat or even decreasing. 

It has been demonstrated by Persic & Salucci (1988, 
1990) and Broeils (1992) that, as the galaxy luminosity 
decreases, the light is progressively unable to trace the ra- 
dial distribution of the gravitating matter (see also Salucci, 
1997). This discrepancy is, in general, interpreted as the 
signature of an invisible mass component (Rubin et.al. 
1980; Bosma 1981). As pointed out by PSS the universal- 
ity of the rotation curves, in combination with the invari- 
ant distribution of the luminous matter, implies an uni- 
versal dark matter distribution with luminosity-dependent 
scaling properties. 

On the theoretical side, recent high-resolution cosmo- 
logical N-body simulations have shown that cold dark 
matter halos achieve a specific equilibrium density profile 
(Navarro, Frenk & White 1996, NFW; Cole & Lacey 1997; 
Fukushige & Makino 1997; Moore et al. 1998; Kravtsov 
et al. 1998). This can be characterized by one free pa- 
rameter, e.g. M200, the dark mass enclosed within the 
radius inside which the average over-density is 200 times 
the critical density of the Universe. In the innermost re- 
gions the dark matter profiles show some scatter around 
an average profile which is characterized by a power-law 
cusp p - r'^, with 7 = 1 - 1.5 (NFW, Moore et al. 1998, 
Bullock et al., 1999). 

Until recently, due to both the limited number of suit- 



able rotation curves and a poor knowledge of the exact 
amount of luminous matter present in the innermost re- 
gions of spirals, it has been difficult to investigate the in- 
ternal structure of dark matter halos. The situation is 
more favorable for (low surface brightness) dwarf galaxies 
which are strongly dark matter dominated even at small 
radii. The kinematics of these systems shows an univer- 
sality of the dark halo density profiles, but, it results in 
disagreement with that predicted by CDM, in particular 
because of the existence of dark halo density cores. (Moore 
1994; Burkert 1995). The origin of these features is not yet 
understood (see e.g. Navarro, Eke & Frenk 1996; Burk- 
ert & Silk 1997, Gelato & Sommer-Larsen 1999), but it 
is likely that it involves more physics than a simple hi- 
erarchical assembly of cold structures. To cope with this 
observational evidence, Burkert (1995) proposed an empir- 
ical profile that successfully fitted the halo rotation curves 
of four dark matter dominated dwarf galaxies 



Pb{r) = 



[r + ro)(r2 



(1) 



where po and are free parameters which represent the 
central dark matter density and the scale radius. This 
sample has been extended, more recently, to 17 dwarf ir- 
regular and low surface brightness galaxies (Kravtsov et 
al. 1998, see however van den Bosch et al. 1999) which 
all are found to confirm equation (1). Adopting spherical 
symmetry, the mass distribution of the Burkert halos is 
given by 



Mb(r) =4Mo{/n(l + ^)-tg-i(^) 



+ \ln 



(2) 



■ 3.2/?^, with the exponential disk length-scale 



with Mg, the dark mass within the core given by: 

Mo = l.&porl (3) 

The halo contribution to the circular velocity is then: 

V^{r)=GMb{r)/r. (4) 

Although the dark matter core parameters ro, po and Mg 
are in principle independent, the observations reveal a 
clear correlation (Burkert 1995): 



Mo = 4.3 X 10^ 



7/3 

\^ kpc ) ^® 



(5) 



which indicates that dark halos represent a 1-parameter 
family which is completely specified, e. g. by the core 
mass. 
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Fig. 1. — Synthetic rotation curves (filled circles with error bars) 
and URC (solid line) with its separate dark/luminous contributions 
(dotted line; disk; dashed line; halo.) See PSS for details. 



The analysis of a recently published large sample of RC's 

(Persic & Salucci, 1995) has provided a suitable framework 
to investigate the dark halo density distribution in spirals. 
The starting points of this study are: a) the mass in spirals 
is distributed according to the Inner Baryon Dominance 
(IBD) regime: there is a characteristic transition radius 

RiBD ^ 2i?d ( 220fcm/J ^^^^ which, at r < Ribd-, the lu- 
minous matter totally accounts for the mass distribution, 
while, for r > Ribd, the DM rapidly becomes the dom- 
inant dynamical component (Salucci and Persic, 1999a,b; 
Salucci ct al, 2000; Ratnam and Salucci, 2000; Borriello 
and Salucci, 2000). Then, although the dark halo might 
extend down to the galaxy center, it is only for r > Ribd 
that it gives a non-negligible contribution to the circular 
velocity, b) The dark matter is distributed in a different 



way with respect to any of the various baryonic compo- 
nents (PSS, Corbclh and Salucci, 2000), and c) The HI 
contribution to the circular velocity, for r < Ropt, is neg- 
ligible (e.g. Rhee, 1996 ; Verheijen, 1997). 

The main aim of this letter is to expand the above re- 
sults to derive the luminosity-averaged density profiles of 
the dark halos and to relate them with the Burkert profiles. 
Section 2 presents the analysis of a homogeneous sample 
of about 1100 rotation curves in which the dark halo con- 
tribution to the circular velocity is first derived and then 
matched to the Burkert halo mass models. In section 3 
we discuss the results. We take Hq = 75km/s/Mpc and 
Oo = 0.3, however no result depends on these choices. 

2. DM HALO PROFILES IN SPIRALS 

PSS (see also Rubin et al. 1982) have derived, from 
15000 velocity measurements of 1000 rotation curves (RC), 
Vsynij^^, j^), the synthetic rotation velocities of spirals 
sorted by luminosity (Figure 1, with Li the I-band lumi- 
nosity and = io-(^'+2i-9)/5). Remarkably, individ- 
MflZRC's have a very small variance with respect to the cor- 
responding synthetic curves (PSS, Rhee 1996, Rhee & van 
Albada 1996, Roscoe 1999, Swaters, 1999): spirals sweep a 
very narrow locus in the RC-profile/amplitude/luminosity 
space. On the other hand, the galaxy kinematical proper- 
ties do significantly change with luminosity (e.g. PSS), 
so it is natural to relate the mass distribution with 
this quantity. The whole set of synthetic RC's has 
been modeled with the Universal Rotation Curve (URC), 
yuRc{f / Ropt^Li / L^) which includes: (a) an exponential 
thin disk term: 



yluRci.x) = 1.28 f3V,lt {loKo - hK^)\i.e. 
and (h) & spherical halo term: 



(1-/3) (l + a^) 



(x2 



(6) 



(7) 



with X = r/Ropt, P = {Vd,uRc{Ropt)/Voptf, Vopt = 
V{Ropt) and a the halo core radius in units of Ropt- At 
high luminosities the contribution from a bulge component 
has been also considered (Salucci and Persic, 1999b). 
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Let us stress that the halo velocity term of eq. (7) docs 
not bias the mass model: it can account for maximum- 
disk, solid-body, no-halo, all-halo, CDM and core-less halo 
mass models. In practice, the values of the free parameters 
a and /? that result from fitting the URC 



URC 



(x) = y^URC (x, P, a) + VluRc (x, 13) (8) 
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to the synthetic curves Vgyn select the actual model. 
Adopting f3 ~ 0.72+ O.Ulog{^) and a ~ 1.5(^)1/^ 
(see PSS) or, equivalently the corresponding 



a{P) 



(i = (3{log Voj>t) 



(9) 



that we have plotted in Fig. (2), the URC mass models 
reproduce the synthetic curves Vsyn{r) within their r.m.s. 
(see Fig. (1)). More in detail, at any luminosity and ra- 
dius, |T^(7flC — Ksynl < 2% and the la fitting uncertainties 
on a and (3 are about 20% (PSS). 

We then compare the dark halo velocities obtained with 
cq. (7) and (9), with the Burkert velocities Vf,(r) of eqs. 
(2)-(4). We leave po and tq as free parameters, i. e. we 
do not impose the relationship of eq. (5). The results are 
shown in Fig (3): at any luminosity, out to the outermost 
radii (~ 6Rd), ^h('') is indistinguishable from Vh,uRc{i^)- 
More specifically, by setting Vh,uRc{r) = Vb{r), we are 
able to reproduce the synthetic rotation curves Vsyn{r) at 
the level of their r.m.s. For r >> 6Rd, i.e. beyond the 
region described by the URC, the two velocity profiles pro- 
gressively differ. 
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Fig. 3. — URC-halo rotation curves {filled circles with error 
bars) and the Burkert model {solid line). The bin magnitudes are 
also indicated. 



The values obtained for tq and po for the URC agrees 
with the extrapolation at high masses of the scaling law 
p oc r^^^'^ (Burkert, 1995) established for objects with 

^The virial halo mass is given by; M200 = 200 X 47r/3pc^J200^o( 
see e.g. Bullock et al., (1999); the critical density is defined as: pc = 



core radii rg ten times smaller (see Fig 4). Let us no- 
tice that the core radii are very large: tq >> Rd so that 
an ever-rising halo RC cannot be excluded by the data. 

Moreover, the disk-mass vs. central halo density relation- 

2/3 

ship Po cx ' , found in dwarf galaxies (Burkert, 1995), 
where the densest halos harbor the least massive disks, 
holds also for disk systems of stellar mass up to IO^^Mq 
(see Fig 4). 

The above relationship show a curvature at the highest 
masses/lowest densities that can be related to the exis- 
tence of an upper limit in the dark halo mass Af2oo ^ which 
is evident by the sudden decline of the baryonic mass func- 
tion of disk galaxies at M^""" = 2 x lO^Mg (Salucci and 
Persic, 1999a), that implies a maximum halo mass of 



^"200 



(10) 



where r^o and f7f, ~ 0.03 (e.g. Buries and Tytler, 1998) 
are the matter and baryonic densities of the Universe in 
units of critical density. From the definition of M200, by 
means of eq. (2) and (3), we can write M200 in terms of 
the "observable" quantity Mq: M200 — V^h- For (^0, z) = 
(0.3,3), ij ~ 12; notice that there is a mild dependences 
of 77 on z and f2o which is irrelevant for the present study. 
Combining eq. (3) and (10) we obtain an upper limit for 
the central density, po < Ix lO^^'^ {ro/kpc)~^g/crrV^, which 
implies a lack of objects with po > 4 x 10~^^.g/cm'^ and 
> SOkpc, as is evident in Figure 3. Turning the ar- 
gument around, the deficit of objects with Md ~ M™°^ 
and Po > 4 X 10~^^ g/cm^, suggests that, at this mass 
scale, the total-to- baryonic mass fraction may approache 
the cosmological value il/fif, ~ 10. 



3. DISCUSSION 

Out to two optical radii, the Burkert density profile re- 
produces, for the whole spiral luminosity sequence, the 
DM halos mass distribution. This density profile, though 
at very large radii coincides with the NFW profile, ap- 
proaches a constant, finite density value at the center, in 
a way consistent with an isothermal distribution. This is 
in contradiction to cosmological models (e.g. Fukushige 
and Makino 1997) which predict that the velocity disper- 
sion fj of the dark matter particles decreases towards the 
center to reach cr ^ for r — » 0. After the result of this 
study, the dark halo inner regions, therefore, cannot be 
considered as kinematically cold structures but rather as 
"warm" regions with size tq oc Pq^'^ ■ The halo core sizes 
are very large: tq 4 — IRd- Then, the boundary of the 
core region is well beyond the region where the stars are 
located and, as in Corbelh and Salucci (2000), even at the 
outermost observed radius there is not the slightest evi- 
dence that dark halos converge to & p (or a steeper) 
regime. 

1 + z''^)g{z) with z the formation redshift, -R200 the virial radius, for g{z) 
■ 3/(87r)G-lH2. 
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Fig. 4. — top Disk mass (in solar units) vs central halo density 
po (in g/cm^) for normal spirals {filled circles). The straight line 
is the extrapolation to high luminosities of the relation of dwarfs. 
bottom Central density vs core radii (in kpc) for normal spirals 
{filled circles), compared with the extrapolation of the relation- 
ship of dwarfs {dotted line, the point with errorbar represents a 
typical object of Burkert, 1995). The solid line is the eye-ball fit; 
po = 5 X 10-'^*r~^^^e-(''o/^''f g/cm^ , The effect of a limiting halo 
mass is also shown {dashed line). 



We find that the dark halos around spirals are essentially 
an one-parameter family. It is relevant that the order pa- 
rameter (the central density or the core radius) correlates 
with the luminous mass (see Fig 4). We do however not 
know how it is related to the global structural properties 
of the dark halo, like the virial radius or the virial mass. 
In fact, the RC out to 6Rd is completely determined by 
the core parameters, i.e. the central core density and the 
core radius, both of which are not defined in the CDM 
scenario. 

The location of spiral galaxies in the parameter space 
of virial mass, halo central density and baryonic mass is 
determined by different processes that occur on different 
scales and at different red-shifts. Yet, this 3D space de- 
generates into a single curve (see Figure 4, remind that: 
Po = 5i Af20o/fo and remind that: Af<j = (3V^ptRopt) 
which describes the dark-luminous coupling. 

Let us discuss the limitations of the present results. 
First, here we have considered the luminosity dependence 
of the dark halo structure. Although this is probably the 
most relevant one, other dependences (Hubble type and 
surface brightness) should also be investigated. Moreover, 
the existence of a (weak) cosmic variance in the halo struc- 
tural properties cannot be excluded until we analyze indi- 
vidual objects (Salucci, 2000, Borriello and Salucci, 2000). 
Secondly, we have derived the profile of DM halos out to 
about six disk-scale lengths, i.e. out to a distance much 
smaller than the virial radius. To assess the global valid- 
ity of the proposed mass model data at larger radii are 
obviously required. 
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